Surface Layering in liquids: The role of surface tension. 
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Recent measurements show that the free surfaces of hquid metals and aUoys are always layered, 
regardless of composition and surface tension, a result supported by three decades of simulations and 
theory. Recent theoretical work claims, however, that at low enough temperatures the free surfaces 
of all liquids should become layered, unless preempted by bulk freezing. Using x-ray reflectivity 
and diffuse scattering measurements we show that there is no observable surface-induced layering 
in water at T=298 K, thus highlighting a fundamental difference between dielectric and metallic 
liquids. The implications of this result for the question in the title are discussed. 

PACS numbers: 68.10.-m, 61.10.-i 
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The free surface of liquid metals and alloys were 
demonstrated experimentally over the last few years to be 
layered, i.e. to exhibit an atomic-scale oscillatory surface- 
normal density profile. ^'^'^'^ This is manifested by the 
appearance of a Bragg-like peak in the x-ray reflectivity 
(XR) curve, R{qz), as shown in Fig. ^for Ga.^ The wave- 
vector transfer position of the peak, qpeak is related to the 
layering period d by qpeak = 2Tr/d. The layered interface 
is in a marked contrast with the theoretical description of 
the liquid- vapor interface of a simple liquid. This theory, 
prevailing for over a century, depicts the density profile 
as a monotonic increase from the low density of the va- 
por and the high density of the bulk liquidi^ This view 
was supported by XR measurements on many non metal- 
lic liquids measured over the last two decades, including 
water^ alkanes;^ and quantum liquids^ which showed no 
Bragg-like peaks. However, as discussed below, the mea- 
surements in all of these studies were restricted to the 
small q range, i.e. g^ <gc vr/atomic size, and would not 
have detected surface layering if it existed. 

Early simulations on non-metallic liquids demonstrate 
that atomic layering is ubiquitous near a hard flat surface 
and this has been observed for liquid Gallium^ On the 
other hand, for the liquid vapor interface it is tempting 
to think that the large surface tension 7 of liquid metals 
such as Hg (7 « 500 mN/m), Ga (7 « 750 mN/m), and 
In (7 « 550 mN/m) might be the explanation for the SL 
observed at their surface. This assumption is partially 
mitigated by the observation of SL at the free surface 
of liquid K where 7 « 100 mN/m only. Here we report 
x-ray scattering results showing that the free surface of 
water, which has nearly the same surface tension as K, 
does not exhibit SL features in the reflectivity proflles, 
thereby suggesting that surface tension by itself does not 
explain SL. This conclusion rests on the validity of the 
capillary wave theory (discussed below) . 

Rice et alM> first predicted SL in liquid metals three 
decades ago. They argued that the layered interface 
structure for liquid metals is a consequence of the strong 
dependence of the effective ion-potential energy on the 



steeply varying electron density across the liquid/vapor 
interface. At the low density vapor phase the electrically 
neutral atoms interact through a van der Waals interac- 
tion only. In the metallic, higher density, liquid phase the 
electrons are delocalized, and the much more complex in- 
teractions involve an interplay between a quantum Fermi 
fluid of free electrons and a classical liquid of charged ion 
cores. Rice concludes that this substantive change in 
the effective-ion-potential stabilizes the short range sur- 
face fluctuations with the result that the atoms near the 
surface form a layered structure J»l Calculations employ- 
ing the glue model of metallic cohesion support these 
conclusions li^ By contrast, Soler et alm^ claim that SL 
is solely due to the formation of a dense layer at the sur- 
face. This layer is not restricted to metallic liquids, but 
may form also in nonmetallic liquids due to nonisotropic 
interactions, such as remnant covalent bonding in liq- 
uid S:^'^ and the highly directional interactions in liquid 
crystalsfi^ or in non-uniform Lennard- Jones fluids with 
unbalanced attractive forces ii^ Additionally, the SL has 
been found in colloidal systemsji^ which suggests that 
layering phenomenon may be more fundamental than 
previously thought. More recently, based on extensive 
simulations, Chacon et al}- argued that SL does not re- 
quire the many-body, delocalize-electron interactions of 
a liquid metal at all. Rather, according to them, SL is 
a universal property of all liquids at low enough tem- 
peratures, T < Tc/ a, whenever not preempted by bulk 
solidification. Here Tc is the critical temperature of the 
liquid, and a « 4 — 5. 

The dichotomy between these two views could be re- 
solved, in principle, by XR measurements on selected 
materialsi^ to see which exhibit, or not, a layering peak 
at some qpeak- Unfortunately, practical considerations 
limit the number of elemental liquids that can be studied 
to a relatively few. One of the major problems is that 
the measurable q^ range is more often than not limited 
to values much less than qpeak by the strong off-specular 
diffuse scattering caused by thermal capillary waves. The 
effect of the capillary waves is to induce a surface rough- 
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FIG. 2: Schematic description of the experimental setup. 



FIG, 1: X-ray reflectivity from free liquid surfaces of the in- 
dicated samples. Points - measured [R{qz)], lines - calculated 
for an ideally flat and steplike interface {Rpi^z)]- The inset 
shows the ratio of the two. 



ness u ~ A/r/7, where 7 is the surface tension. The 
consequence of this, which is shown in Fig, ^ for three 
hquids at room temperature, is to reduce the reflectiv- 
ity R{qz) below that of the theoretical Fresnel reflectivity 
from an interface with an idealized flat, steplike surface- 
normal density profile. For low-7 liquids such as water (^ 
70 mN/m) and K (~100 mN/m) the reduction is signifi- 
cant. For Ga, where 7 k, 750 mN/m the effect is almost 
negligible at room temperature. Even that small reduc- 
tion is almost completely offset at room temperature by 
the SL effect that peaks at q^ ~ 2.5 A^^. At higher 
temperatures, however, the effect is quite prominent.^ 

Although the rapid falloff in R/Rp of water prevents 
its measurements out to Qz ~ Qpeak = 2,0 A^^, our recent 
studies of liquid K^a^ demonstrated that if surface layer- 
ing is present its signature can still be observed clearly 
at g^ << qpeak even for low 7 < 100 mN/m liquids. This 
is accomplished by carefully accounting for the effects 
of capillary waves, based on diffuse x-ray scattering (DS) 
measurements. We present here an x-ray study of the sur- 
face structure of water over the most extended gz-range 
published to date, and including diffuse scatteringL^S. 
For the present measurements the intrinsic surface struc- 
ture factor of water can be extracted directly from the 
raw R{qz) without resorting to any structural model for 
the interface. Comparison between the water surface 
structure factor, for which there is no evidence of SL, and 
that of K and Ga suggests that surface tension is not the 
dominant cause of the SL observed in liquid metals. 

X-ray measurements were carried out on the CMC- 
CAT liquid surface diffractometer, APS, Argonne Na- 
tional Laboratory, at a wavelength of A = 1.531 A. 
The purified water sample was contained in a Langmuir 
trough-^ mounted on the diffractometer. The surface 
was periodically swept with a teflon barrier, monitoring 



7 with a film balance, to ensure a clean surfaced 

Both XR^^ and DSSi are well documented techniques. 
Experimental geometry setup employed in the experi- 
ments described in this work is shown in Fig. 12 For x- 
rays striking the surface at a grazing angle a and detected 
at an output angle (3 in the specular plane, the surface- 
normal (z) and in-plane surface-parallel {y) components 
of the wavevector transfer are qz = {2tt/X) (sin a + sin /?) 
and qy — (27r/A) (cos a — cos/3). In a XR measurement 
both angles are varied, keeping a — P, so that qy — 0, 
and the reflected signal, measured vs qz and divided by 
the incident intensity /q, yields R{qz). In a DS measure- 
ment, (3 is varied for a flxed a, and hence both qz and qy 
vary. 

Within the Born 

R{qz)/RF[qz) = Mqz)?W{i^,qz), 
(Poo) V[d(p(^))/d2]exp(zqzz)dz is 
structure factor of the liquid-vapor 
is the intrinsic (i.e. in the absence of capillary wave 
smearing) surface-parallel-averaged electron density pro- 
flle along the surface-normal z direction, and W{ri,qz) 
accounts for the smearing of the intrinsic density profile 
by capillary waves. The aim of our XR measurement is 
to determine ^{qz) to observe a possible layering peak. 
Since both '^{qz) and W{ri,qz) depend on qz, extracting 
^{qz) directly from the measured R{qz) is possible as 
|$(gz)P = R[qz)/W{-q,qz) only if W{r],qz) is known 
independently. For this we use the DS data. 

The measured DS, shown in Fig. |31 is given by theory 

as Ids = J['i<^{qx-q'x,qv-qy,qz-q'z)/<i^]'i^{q'xiq'y,q'z), 

where duj is the angular resolution of the diffractometer 
and 



approximation^ 
where ^{qz) = 
the conventional 
interface, {p{z)) 
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is the scattering cross-section. Here qc is the critical angle 
for total external reflection of x rays, qmax ~ t^/S, is the 
upper cutoff for capillary wave contributions, with ^ of 
order of the atomic diameter, and 77 = {ksT) / {2'K^)q^. 
Two subtle complications are encountered in analyzing 
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FIG. 3: Comparison of measured diffuse scattering with cap- 
illary wave theory predictions for the angles of incidence a 
(top to bottom): 2.1°, 2.8°, 3.5°, 4.2°, 5.0°, 5.7°, 6.0°, and 
6.4°. The last data set corresponding to 7.1° shown in the 
inset no longer exhibits a distinguishable specular peak. The 
qz values corresponding to the specular condition qy — Q A~^ 
are 0.3, 0.4, 0.5, 0.7, 0.8, 0.85, 0.9, and 1.0 A"\ respectively. 
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FIG. 4: Comparison of the structure factor squared | $(gz) P 
for water (crosses), liquid potassium (squares), and liquid gal- 
lium (triangles). The wavevector q^ is normalized to the ex- 
pected position of the layering peak q^eak of each sample. The 
inset shows the data on an extended scale. For discussion see 
text. 



the DS data. First, nonsurface DS contributions (scatter- 
ing from the bulk, sample chamber windows, etc.) can 
significantly distort the shape of the DS scans. These 
background contributions are measured by offsetting the 
detector by A8 = ±0.3° from the plane of incidence, and 
are already subtracted from the data shown in Fig. 13 
Second, in the fixed-a DS scans, carried out by scan- 
ning /?, 92, and thus |$((72)p, also vary. Fortunately, 
5qz ~ Sqy/P and for the range of qy displayed in Fig. O 
and typical qz values, the changes are small enough to 
be neglected, e.g. a maximal Sqz ~ 0.06 A at the 
largest /3 for a = 3.5°. Thus, it is a good approxi- 
mation for each of the DS scans to treat $((?z) as a 
fixed function of a. | $((7z) p is then obtained by di- 
viding the measured DS and XR curves by W{ri,qz) — 
J[Ao/{87TsmaMRFiqz){v/ql-'')q;2.doj^ 

The theoretical curves calculated using Eq. Q] with 
T ~ 298 K and 7 = 72 mN/m are shown as lines in Fig. 
13 As qz increases, so do both f] oc qz and the intensity of 
the off-specular power-law wings relative to that of the 
specular peak at qy = 0. The curve at qz = 1 A~^ demon- 
strates the capillary-wave-imposed limit where the spec- 
ular signal at qy = 0, which contains the surface structure 
information, becomes indistinguishable from the DS sig- 
nal at ±qy > 0. In principle, this limit arises from the 
fact that for 77 > 2 the singularity at gj^ = in da/dfl 
vanishes and there is no longer any criterion by which 
the surface scattering can be differentiated from other 
sources of diffuse scattering. In practice, the fact that 
the projection of the resolution function on the horizon- 
tal x-y plane is very much wider transverse to the plane 
of incidence than within the plane of incidence reduces 



this limit to a value closer to ry ~ 1.^4 Figure 01 exhibits 
excellent agreement between the theoretical DS curves 
calculated from Eq. ^with the measured DS over several 
decades in intensity and one decade in 77, without any 
adjustable parameters. This confirms the applicability of 
the capillary wave theory for the surface of water over 
the qz range studied here, < gz < 0.9 A~^. Measure- 
ments of diffuse scattering for small qz (i.e. small 77) to 
values of the surface parallel component of the wave vec- 
tor transfer q^ of the order of 7r/atomic size have been 



done by Daillant et. 
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by moving the detector out 



of the plane of incidence for grazing incident angles. At 
such large wave vectors the observed scattering must be 
interpreted as the superposition of scattering due to sur- 
face capillary waves and bulk diffuse scattering from the 
liquid below the surface. The separation of these two 
contributions is rather subtle and was only accomplished 
through nontrivial calculations of the noncapillary terms. 
Eventually Daillant et. al concluded that as the value of 
qx approaches the atomic scale (i.e., 7r/atomic size) the 
surface tension varies with qx ^ In principle, this disper- 
sion should affect the sum rule that is the origin of the 
l/gjr, in the differential cross section in Eq.(^. On the 
other hand, considering that the dispersion is a relatively 
small effect that will only change the argument of a log- 
arithmic term, its effect on the present analysis can be 
neglected. 

I'&l'Zz)!^ is then obtained directly from the measured 
XR curve, as R{qz)/W{r],qz) as discussed above. It is 
shown in Fig. ^ (circles) along with previously measured 
results for K (squares) and Ga (triangles). 

The rise of the Ga |$(gz)P to '^100 at qpeak due to 



layering can be clearly seen on the inset of Fig. 01 For 
K, the capillary- wave-imposed limit only allows obtain- 
ing |<I>(g2)|^ for Qz < O.Sqpeak- Nevertheless, the value of 
$(gz)| for K starts to deviate from unity for values of 
Qz/lpeak ~ 0.3. Furthermore, over the range for which it 
can be measured it is basically identical to the structure 
factor of Ga. This is a clear indication that the surface of 
liquid K has essentially the same SL as that of Ga, which 
is also nearly identical to that of the other liquid met- 
als that have been studied to date, e.g., InSl and SnSSi. 
For water, however, no deviation of |<i>(q2)p from unity is 
observed even at the highest measurable Qz/qpeak ~ 0.5. 
This suggests that surface-induced layering does not oc- 
cur at the surface of water. The different behavior, in 
spite of the similar 7 of water and K, leads to the con- 
clusion that the surface layering in K, and by implication 
in other liquid metals, is not merely a consequence of its 
surface tension. 

The absence of layering in water, and its presence in 
potassium, seems at first sight to corroborate the claim 
of Rice et alw^ that layering is a property arising from 
the metallic interaction of the liquid. On the other hand, 
Chacon et ali^, who maintain that surface-induce layer- 
ing is a general property of all liquids, regardless of their 
interactions, predict that layering should occur only at 
temperatures T/Tc < 0.2, where Tc is the critical temper- 
ature of the liquid. Although supercooling is often possi- 
ble, the practical limit for most reflectivity measurements 
is the melting temperature T^. Thus the smallest T/Tc 
for any liquid is on the order of T^/Tc. For liquid metals 
Tm/Tc « 0.15 (K), 0.13 (Hg), 0.07 (In), 0.066 (Sn), and 
0.043 (Ga). Since these values are < 0.2, by Chacon's 
criteria surface layering is expected, and indeed demon- 
strated experimentally to occur, in all of thempiiSii2iSZiS& 
By contrast, for water, where Tm/Tc = 0.42 ^ 0.2, 
Chacon's criteria predict that the appearance of surface 
layering is preempted by bulk freezing, and thus no SL 
should occur at room temperature, as indeed found here. 
In summary, we have shown that the surface of water 
does not exhibit SL even though its surface tension is 
not significantly different from that of liquid K for which 
SL was observed. Thus the atomic layering of liquid K 
cannot be the result of surface tension alone. Although 



this seems to support Rice's argument that the metallic 
phase is essential for SL, we note that SL is exhibited by 
both liquid crystalsi^ and other intermediate-size organic 
molecules. ^^ Consequently, there can be other criteria for 
SL, in addition to those proposed by Rice. One possibil- 
ity, proposed by Chacon et. al, is that SL should be 
ubiquitous for all liquids that can be cooled to tempera- 
tures of the order of 0.2rc. Unfortunately, this is difficult 
to explore experimentally since suitable liquids are rather 
scarce. For example, liquid noble gases, the archetypical 
van der Waals liquids, have all T^jT^ > 0.55. For liquid 
helium the ratio is lower but in view of both the pres- 
ence of the superfluid transition and the low scattering 
cross section this is a difficult system to study. X-ray 
reflectivity measurements at 1 K, i.e., Tm/Tc ~ 0.2, did 
not exhibit evidence for SL^. Similarly, the polar liq- 
uids oxygen and fluorine have T„i/Tc of only 0.38 and 
0.37. There are, however, some low melting organic 
van der Waals liquids, e.g. propane and 1-butene, that 
have 0.2 < T„i/Tc < 0.25, and may be suitable for ad- 
dressing this issue. It would be also very interesting to 
measure the surface structure of materials such as Se 
{T,n/Tc = 0.28), Te (0.3), Sc (0.28), and Cd (0.22), for 
which Tjn/Tc > 0.2. If it were possible to probe surface 
layering in these high Tm/Tc materials such experiments 
could test Chacon's argument. Unfortunately, these ma- 
terials have relatively high vapor pressure and that pre- 
cludes the use of UHV methods for insuring atomically 
clean surfaces. The surfaces of some alkalis and mercury 
can be kept clean without UHV methods. For most of 
the experimentally accessible metals this is not true. 
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